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impurities 

A. P. Misra^'Eland N. C. Adhikary2'E3 

Department of Mathematics, Siksha Bhavana, Visva-Bharati University, Santiniketan-731 235, 
India. 

^■^ Material Sciences Division, Institute of Advanced Study in Science and Technology, Paschim Boragaon, 
Garchuk-781 035, Guwahati, Assam, India 

(Dated: 11 November 2011) 

The nonlinear propagation of large amplitude dust ion-acoustic (DIA) solitary waves (SWs) in an ion-beam 
plasma with stationary charged dusts is investigated. For typical plasma parameters relevant for experiments 
[J. Plasma Phys. 60, 69 (1998)], when the beam speed is larger than the DIA speed [v^q > l.lcg), three stable 
waves, namely the "fast" and "slow" ion-beam modes and the plasma DIA wave are shown to exist. These 
modes can propagate as SWs in the beam plasmas. However, in the other regime (cs < Vbo < l-Tcg), one of 
the beam modes when coupled to the DIA mode may become unstable. The SWs with positive (negative) 
potential may exist when the difference of the nonlinear wave speed (M) and the beam speed is such that 
1.2 < M — Vbo ^ 1.6 (Af — > 1.6). Furthermore, for real density perturbations, the wave potential (> 0) is 
found to be limited by a critical value which typically depends on M, Vbo as well as the ion/beam temperature. 
The conditions for the existence of DIA solitons are obtained and their properties are analyzed numerically 
in terms of the system parameters. While the system supports both the compressive and rarefactive large 
amplitude SWs, the small amplitude solitons exist only of the compressive type. The theoretical results may 
be useful for observation of soliton excitations in laboratory ion-bcam driven plasmas as well as in space 
plasmas where the charged dusts play as impurities. 



PACS numbers: 52.35.-g, 52.30.Ex, 52.25.Vy 



I. INTRODUCTION 

Large amplitude solitary waves (SWs) in plasmas with 
a high energy ion beam have been observed at vari- 
ous space plasma environments, e.g., the Earth's mag- 
netopause, in the Van Allen radiation belts as well as 
in the auroral zoneii^. These nonlinear waves may be 
driven by the momentum exchange between the elec- 
trons and protons, or between different ion populations in 
multi-ion plasmas. Propagation of such ion-acoustic (lA) 
SWs (lASWs) in magnetized or unmagnetized collision- 
less plasmas under different physical situations has been 
of considerable interest in recent years^"— . Using different 
methods, many authors have studied the behavior and 
characteristics of lA solitons both theoretically^i^i^"— 
and experimentally^ii^^— in ion-beam plasmas. 

It has been found that the presence of energetic 
charged particles like ion beam in plasmas can signifi- 
cantly modify the propagation behaviors of SWsi^. The 
latter with negative potentials have been found in the 
vicinity of ion beam regions of the auroral zone in the 
upper atmospherei^. The spacecraft observations in the 
Earth's plasma sheet boundary indicate that the elec- 
tron and ion beams can also drive the broadband electro- 
static waves^. Furthermore, the low-temperature plas- 
mas containing charged dusts are ubiquitous in various 
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space plasmasi^, laboratory dcvicc3.i2, as well as in indus- 
trial processes^. The presence of such negatively charged 
dusts can significantly influence the collective behaviors 
of plasmas^ii^. 

On the other hand, it has been shown that when 
an ion beam is injected into an unmagnetized plasma, 
three stable normal modes, namely the "fast" and 
"slow" ion-beam modes and the plasma lA wave exist 
in plasmasi^iii. Furthermore, depending on whether the 
beam speed is of the order of or larger (e.g. 2 times) 
than the lA speed, the lA waves when coupled to slow 
beam mode may become unstable. Such ion-ion insta- 
bility has been confirmed experimentally by Gresillon et 
aJi^. Thus, there may exist either three stable modes 
or one unstable and two stable modes. The nonlinear 
wave evolution of such three modes was investigated by 
Yajima et al^ in an ion-beam plasma. They showed 
that each of these modes can propagate as small ampli- 
tude Korteweg-de Vries (KdV) solitons when the beam 
density is smaller than the electron density and the am- 
plitude is smaller than a critical value. In an another 
work, Nakamura et al^ had experimentally shown that 
the propagation of large amplitude lASWs (compressive 
type) is possible for beam speed larger than the lA speed 
in an ion-bcam plasma with two groups (high and low- 
temperature) of electrons. However, it seems that the 
effects of adiabatic positive ions and ion beam as well 
as their finite temperatures on the propagation of arbi- 
trary amplitude dust lASWs (DIASWs) in a plasma with 
charged dust impurities, have not yet been considered in 
detail. 
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In this paper, we study the properties of three hnear 
eigenmodes, which may propagate as stable DIASWs. 
For typical plasma parameters relevant for experimen- 
tal conditionsi^, one of the modes may become unstable 
when the beam speed to dust ion-acoustic (DIA) speed 
ratio lies in 1 < Vbo/cg < 1.7. The conditions for the 
existence of large amplitude DIASWs are analyzed by 
pseudopotential approach and their properties are stud- 
ied numerically for fast modes. Depending on the pa- 
rameters, we show the existence of compressive as well 
rarefactive DIASWs. The latter exist when the ion to 
electron temperature ratio is of the order ^0.1 or more. 
On the other hand, the DIASWs with small amplitudes 
are shown to propagate only of the compressive type. 



II. BASIC EQUATIONS AND THE DISPERSION 
RELATION 

We consider an unmagnctized collisionless plasma com- 
posed of adiabatic positive ions, positive beam ions, 
Boltzmann distributed electrons and negatively charged 
(stationary) dust grains in the background plasma. The 
finite temperatures of both the beam and plasma ions 
together with an equilibrium flow of the beam are con- 
sidered. The normalized equations read 
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where rij, Vj and Pj respectively denote the number 
density, speed and the pressure of singly charged beam 
(j = b) and plasma ions (j = i), normalized by the 
unperturbed number density rijo, the ion- acoustic speed 
Cs = ^JkBTe/rrii and njoksTj. Here kB is the Boltzmann 
constant, Tj is the temperature for j-th species (j = e 
for electrons) and mj is the particle mass with nii ~ mi,. 
Furthermore, is the DIA wave potential normalized by 
ksTe/e, where e is the elementary charge. The space and 
time variables are respectively normalized by the Debye 
length A_D = ksTe / Anniae^ and the ion plasma period 
■ ~^ — 1/ \/ 47rriioe^ / . The overall charge neutrality 



pi 



condition in the background plasma reads 



/i5 = 1 + /id, 



(5) 



where = n^o /^^eo and /id = ZdUdo /neo in which Udo is 
the dust number density and Zd is the number of elec- 
trons on the dust grain. In Eq. ([2]), aj = Tj/Tg is the 
temperature ratio for j = i, h. We have neglected the 
electron inertia, since the electron thermal speed is much 
larger than the ion/beam speed. 



It has been theoretically and experimentally shown 
that when beam ions are injected into an unmagnctized 
plasma, three longitudinal modes involving ion motions, 
namely an ion-acoustic wave and the fast and the slow 
space charge waves in the beam, can propagatei^ill. In 
order to identify those modes in our dusty beam plasma 
system we linearize the basic equations and assume that 
the perturbations vary as ^ exp(ifcx — iut), where ui is 
the angular frequency normalized by uipi and k is the 
wave number normalized by A^^. Thus, we obtain the 
following dispersion relation. 



(6) 



where = 1 + /c^ and wto is the speed (drift) of beam 
ions in equilibrium. In absence of the ion beam, DIA 
mode propagates with the phase speed given by 



k 



(7) 



However, in absence of plasma ions the fast (F) and slow 
(S) modes propagate in the beam with the phase speeds. 



k 



VbO ± A/ScTb 



fj-b 



(8) 



When the beam speed is much larger than the ion- 
acoustic speed, the beam and background ions are not 
strongly coupled. The background ions support slightly 
modified ion-acoustic waves and the beam ions support 
ion-acoustic waves (F and S modes) whose phase speeds 
arc shifted by the beam speed Vbo- Equation ([8|) shows 
that the phase speed of the fast (slow) modes is larger 
(smaller) than the equilibrium beam speed VbQ. Now, the 
dispersion equation © can be expressed as a polynomial 
equation in uj of degree 4, and when the phase speed is 
much larger than the DIA speed, the dispersion relation 
gives real roots of the wave frequency. However, in the 
opposite case, the coupled wave modes can be stable for 
VbO > Vbc where 



"be 



and unstable when 



3(cti + ab) + 



X2 ' 



SfJb 



3<7iJIb 



< vIq < 3{ai + ab) + 



(9) 



,(10) 



For typical laboratory plasma parameters'^ = 0.5, 
CTfc = 0.01, /Xi -I- //f, = 1.4 and k = 0.1, the critical value 
of VbQ above which the modes are stable is Vbc ~ 1.7, and 
for unstable modes we have 1 < Vbo < 1.7. Thus, it is 
reasonable to consider higher values of the beam speed 
(i.e. roughly greater than 2 times the DIA speed) in or- 
der to avoid ion-ion instability^. The latter has been 
confirmed in experiments'^. 

Next, we numerically examine the dispersion relation 
^ for the three wave modes which propagate along the 
beam direction. The results are shown in Fig. 1 where a 
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FIG. 1. (Color online) The dispersion equation ([6]) is contour plotted against the normalized wave number (k) and the 
frequency (cj) for fast (F), slow (S) and DIA modes. The parameter values are (a) i}bo = 5, ai = 0.05, at = 0.01, fid = 0.4 with 
different /i^: /li — 0.2 [thin (red) line] and fn = 0.7 [thick (blue) line]; (b) vto = 4, /^i = 0.2, at = 0.01, Hd = 0.4 with different 
ai: ai = 0.05 [thin (red) line] and ai = 0.5 [thick (blue) line]; (c) Hi — 0.2, ai = 0.05, ai, = 0.01, fid = 0.4 with different vi,o: 
Vbo = 4 [thin (red) line] and vto = 5 [thick (blue) line] and vto = 4, ai = 0.05, at — 0.01, fii — 0.2 with different ftd'. fid ~ 0.4 
[thin (red) line] and fid ~ 0.8 [thick (blue) line]. 



quartic equation in uj is contour plotted against the wave 
number. The behaviors of the modes are found similar 
as experimentally observed modes but in different beam 
plasmasi^. It clearly shows how the frequencies of the 
DIA wave as well as the F and S modes get modified 
by the increase of the beam or ion density, the beam 
speed, the ion temperature as well as the percentage of 
impurity in the background plasma. We find that as /i; 
increases [Fig. 1(a)] the phase speed of the F-modes (S- 
modes) decreases (increases) and that for the lA mode 
increases. However, the ion-temperature only modifies 
the DIA mode as in Fig. 1(b). Figure 1(c) shows that the 
effect of the beam speed is to enhance the phase speeds 
of both the F and S modes. In contrast to Fig. 1(a), 
the enhancement of the charged dust concentration fid 
results into the increase (decrease) of the phase speed of 
F (S) modes, whereas the DIA modes remain almost un- 
changed by the effect of fid- In the following section we 
determine the conditions for which these stable modes 
would indeed propagate as SWs, and analyze their prop- 
erties with different plasma parameters using the pseu- 
dopotential approach. 



III. LARGE AMPLITUDE SOLITONS: 
PSEUDOPOTENTIAL APPROACH 



Assuming that the perturbations vary in the moving 
frame of reference ^ — x ~ Mt, where M is the nonlinear 
wave speed normalized by (If M be normalized by the 
phase speed of the ion- acoustic waves, it would then be 
called as the Mach number), we obtain from Eqs. (Il])-(l4]) 
the following relation for the densities. 



a/2ct ,0 



$1 



1/2 



(11) 



where 



0, $, 



1 - 20/ (M - vjo) 



2 2 



1 



<7^Q and 



2cr 



/4- 



In our nonlinear theory we will consider M to be larger 
than the beam speed in order to examine whether the 
fast modes propagate as solitary waves^^. Similar anal- 
ysis can also be done for the slow modes with M < Vbo- 
In obtaining Eq. (jlip we have used the boundary condi- 
tions, namely (/)—?► 0, Vi^t (0,Wbo), n-j — >■ 1 and Pj — > 1 
as ^ — >■ ±oo. Furthermore, inspecting on Eq. (jll[) . we 
find that the number densities are real when there exists 
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a critical value 4>c of such that Q < cj) < (pc where 

(^c = min 1^ (1 - a^f , ^{M - vtof (1 - abo f 

(12) 

Note that the above restriction is valid only for large 
amplitude waves with positive potential. However, for 
waves with negative potential, the values of (j) are not 
limited by the ratios, namely fii, fib or Ci, at- Typically, 
for M > Vbo and at < at < I, (f>c assumes the second 
term in the curly brackets in Eq. p2)) . However, this 
value must be considered together with the conditions for 
the existence of SWs. Introducing now the relation PT|) 
into Eq. and integrating it we obtain the following 
energy balance equation for an oscillating particle of unit 
mass at the pseudoposition (j) and pseudotime ^. 



1 
2 



Vi(f>) = 0, 



(13) 



where the pscudopotcntial V is given by 
y(0) = l-exp(0)-Aid0-^ 

3/2 



1 - a' 



3/2 



3a: 



-1/2 



- l-Jl-a 



-1/2 



(14) 



Here /3j = iij{M - Vjo)'^<J^i/2\/2ajo. Equations (HH) and 
p4p are valid for arbitrary amplitude stationary pertur- 
bations like SWs and/or double layers. The conditions 
for the existence of such perturbations can be obtained 
as follows: 

(i) V{Q) ~ 0. This has already been satisfied in obtain- 
ing Eq. (|14p . and by using the charge neutrality condition 
([S]); one can easily verify that dV {(p) / d(t) = at = 0; 

(ii) d'^V{(j))/d(j)'^ < at = 0. This is satisfied when 
the following inequality holds. 



fJ'b 



M2 - 3a, {M - Vbo)^ - 3ab 



< 1. 



(15) 



(iii) V^((/)m 7^ 0) = and dV {(j)jn) / d<j) ^ according to 
whether the SWs are compressive (with positive poten- 
tial, i.e. (f> > 0) or rarefactive (with negative potential, 
i.e. (f) < 0). Here 0,„ represents the amplitude of the 
solitary waves or double layers, if exist. 

In order that the three modes propagate as SWs the 
conditions (i)-(iii) must be satisfied. However, for the 
double layers to exist there must be an additional con- 
dition, i.e., dV{(f>m)/d4> = to be satisfied along with 
(i)-(iii). We numerically analyze the conditions (ii) and 
(iii) for some typical plasma parameters as relevant for 
experiments^. The condition (ii), which results to Eq. 

is presented in the fii — M plane as shown in Fig. 
2 for different sets of parameters. The white (gray or 



shaded) regions show the parameter space where the 
inequality ([T5|) is satisfied (not satisfied). The upper 
boundary curves that separate the white and gray re- 
gions represent the minimum values of Af (> Vbo) for the 
existence of solitary waves or double layers. The upper 
limits of M (with M > Vbo and ab < (Ji < 1) for the ex- 
istence of positive SWs can be obtained by substituting 
= 0, = i (M - Vbof (1 - <Jbof into 1/(0) = [cf. Eq. 



(IT4l) ]. These limiting values of M are represented by the 
thick (blue) curves in Fig. 2. Typically, for M > Vbo, 
Fig. 2 shows that the lower and upper limits of M in- 
crease with the beam speed as well as the charged dust 
concentration. However, the case of ion temperature ef- 
fect [cf. Figs. 2(a) and (c)] is different. Here the range 
of M (> Vbo) remains almost the same with the increase 
of (Ti. Thus, for M > Vbo, SWs may exist in the upper 
parts of the white regions. 

Figure 3 shows the contour lines otV{(j)) = [thin (red) 
curves] and dV[(f))/d(j) = [thick (blue) curves] for dif- 
ferent values of M . There are, in general, wide ranges of 
values of /i; and 0m , for which V = and dV/d4> = can 
be satisfied. However, all the numerical values of = 
are not admissible as values of = 0„i is limited by (f>c 
such that < (f>m < 4>c [cf. Eq. P^ ]. Such (pc changes 
with M or some other parameters ai, ab- Thus, for each 
set of parameters we can obtain a value of (pc below (for 
positive (p) which SWs may exist. Figures 3(a) and (d) 
explain the variation of the charged dust impurity and 
Figs. 3(b) and (c) that due to ion temperature in dif- 
ferent (higher and lower) regimes of (f)m- We find that 
the enhancement of the ion temperature (0.1 < cr; < 0.5) 
can lead to the existence of negative SWs [Figs. 3(b) and 
(c)]. The common parameter regimes (left to or above 
the curves) in Figures 3(a) and (d) where both the con- 
ditions (j> > and dV/dcj) > arc satisfied gives the 
region of the existence of positive SWs, whereas the con- 
ditions <j> < and dV/dcf) < arc satisfied only to the 
left to or down the curves in Figs. 3(b) and (c), imply- 
ing the existence of negative SWs. Note, however, that 
with similar parameter values as in Fig. 3(b) or (c), the 
formation of negative SWs was not reported in beam- 
plasma experiments^^. Since there is no common point 
of intersection of 1/(0) = = dV{4>)/d(j> except at = 0, 
we conclude that the double layer formation may not be 
possible in our system. 

Next, the pseudopotential 1^(0) can be plotted against 
for a given set of physical parameters. We see from 
Fig. 4 that for both positive and negative values of 0, 
1/(0) crosses the 0-axis for the fast modes, implying the 
existence of both positive and negative SWs. These dif- 
ferent crossing points for certain values of the parameters 
give the height or amplitude \4>m\ of the SWs. However, 
the width of the SWs with speed M can be obtained 
from the shape of 0(^) by the numerical solution of Eq. 
p3)) . We can also verify the wave amplitude and width 
[|0m|/|14nin|, whcrc |l/min| is the minimum value of 1/(0)] 
obtained from the plots of 1/(0) with those from the nu- 
merical solution of Eq. ([T^]) . Since the values of are 
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FIG. 2. The parametric regions (white or blank) are shown for d^V/dcp'^ < at = in the {M — fn) planes for (a) vto ~ 4, 
fid ~ 0.4, ai = 0.05 and = 0.01; (b) the same as in (a) but Vbo ~ 5; (c) the same as in (a) but ai — 0.5; and (d) the same as 
in (a) but fj,d = 0.8. The shaded (gray) regions correspond to d^V/dcj)^ > at </) = 0. Thus, for M > vto, SWs may exist only 
in the upper parts of the white regions. The thick (blue) curves represent the upper limits of M for the existence of positive 
SWs. 



limited by (j)c, i.e., < (j) < (pc , and (pc changes with 
M and other parameters, V{4)) can cross the t^-axis only 
for a set of parameters. Wc find that the amplitude of 
the wave increases with the enhancement of the ion con- 
centration until /ii < 0.6 [Fig. 4(a)] as well as the wave 
speed M < 5.4 [Fig. 4(b)], and it decreases with the 
increase of the dust impurity, fid > 0.4 [Fig. 4(d)]. On 
the other hand. Fig. 4(c) shows the existence of negative 
SWs with large negative values of (f>. The amplitudes (in 
absolute value) of these SWs increase with the nonlin- 
ear wave speed M. For the parameters as in Fig. 4(c) 
we find that in contrast to positive SWs in which upper 
limits of M exist, there exist a lower limit of AI (> 5.6) 
above which the SWs with negative potential exist. 

To verify the results as in Fig. 4 (especially the am- 
pHtude and width of the SWs) and to obtain the shape 
of (/){£,), wc numerically integrate Eq. ([T5|) . The results 
are presented in Figs. 5 and 6 for the positive and neg- 
ative SWs respectively. Wc find that for compressive 
SWs (Fig. 5), the amplitude increases and the width de- 
creases with increasing values of M and fit , whereas both 
the amplitude and the width are enhanced by the effect 
of charged dust impurity /x^. Figure 6 shows that the 
effects of M and ai respectively are to increase and de- 
crease the amplitude as well as the width of the negative 
SWs. 



IV. SUMMARY AND CONCLUSION 

The nonlinear evolution of large amplitude SWs in 
a dust contaminated ion-beam driven plasma is inves- 
tigated by using a pseudopotential approach. The dis- 
persion properties of three linear modes, namely the fast 
and slow ion-beam modes and the DIA wave are analyzed 
numerically. The conditions for which these three modes 
propagate as (large amplitude) SWs and their (soliton) 
properties arc studied numerically in terms of the system 
parameters. While the system supports both compressive 
(positive) and rarefactive (negative) large amplitude DT 
ASWs, the small amplitude DIA solitons exist only of the 
compressive type (see Fig. 5). The presence of charged 
dusts significantly alters the existence regions as well as 
the properties of sohtons. We also show that for typical 
plasma parameters as in experiments^^ the formation of 
double layers is not possible. Our results can be summa- 
rized as follows: 

(i) When the beam speed is larger than (~ 2 times) 
the DIA speed, each of the stable modes, namely the 
fast and slow beam mode as well as the DIA mode can 
propagate as SWs in ion-beam plasmas. Otherwise, there 
may exist one unstable beam mode coupled to the DIA 
wave and two stable modes in the plasma. With the en- 
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FIG. 3. (Color online) The contour lines are shown along which V{(j}m) = (thin lines) and dV {(j)m) / d(j) ~ (thick lines) for 
different M as in the figure. The other parameter values for the subplots (a), [(b) and (c)] and (d) are the same as in Fig. 2(a), 
Fig. 2(c) and Fig. 2(d) respectively. The values of (ji = (f)m are such that < 0m < (j>c for positive SWs, where 0c is given 
by Eq. (|12p . The subplots (b) and (c) are shown in two different regimes of where the formation of rarefactive (negative) 
SWs may be possible. The plots (a) and (d) show the regions for the existence of compressive SWs. In all the plots, there is 
no common point of intersection of V{(f)m) = and dV {(j)m) / d(j) = {(j>m 7^ 0) implying that the formation of double layers is 
not possible. The common regions where (j)m > and dV{4>m)/d4> > [left to or above the curves in plots (a) and (d)] are for 
positive SWs and those where (j>m < and dV{4)m)/d4) < [left to or below the curves in plots (b) and (c)] are for negative 
SWs. 



hancement of fXi, while the phase speed of the F- modes 
(S-modes) decreases (increases) that for the DIA mode 
increases. These behaviors arc in contrast to the effect 
of charged dusts fid in which the DIA modes remain al- 
most unchanged. The effect of the beam speed is to en- 
hance the phase speed of both the F and S-modes. In 
the limit of the phase speed larger than the DIA speed, 
the modes propagate as stable waves. However, in the 
opposite limit, the DIA mode, coupled to a beam mode, 
may become unstable in 1 < Vbo < 1.7 for typical plasma 
parameters cr^ = 0.5, Ub = 0.01, /i^ + /i^ = 1.4 at fc = 0.1, 
as relevant for experimental. Thus, it is reasonable to 
consider higher values of Vbo (> 2) for the DIASWs to 
exist. 

(ii) For the perturbations to be real, the wave potential 
(j) (> 0) is to be less than a critical value, which typically 
depends on the nonlinear wave speed, the beam speed as 
well as the beam or the ion temperature. 

(iii) Both the compressive and rarefactive large ampli- 
tude DIASWs may coexist, whereas the small amplitude 
soliton exists only of the compressive type. In order that 
the the SWs with > (0 < 0) may exist, the regime of 



the nonlinear wave speed (M) and the beam speed {vbo) 
is such that 1.2 < M - Vbo < 1-6 (M - vto > 1-6). 

(iv) For large amplitude positive SWs, the effects of 
M and fii arc to enhance the wave amplitude and to 
reduce the width. However, both the amplitude and the 
width may be increased (in magnitude) by the charged 
dust concentration fid- These behaviors are similar to 
the case of large amplitude negative SWs by the effect of 
M. 

The theoretical results could be useful for soliton ex- 
citation in laboratory ion-bcam driven plasmas as well 
as in space plasmas where negatively charged dusts are 
considered as impurities. 
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FIG. 4. (Color online) The pseudopotential V{(f>) is plotted against (j) to show the existence of positive [subplots (a), (b) and 
(d)] and negative [subplot (c)] SWs. Subplot (a): The solid, dashed, dotted and dash-dotted lines respectively correspond to 
fii = 0.2, 0.4, 0.6 and 0.7 for a fixed A'l — 5.2. The other parameters are as in Fig. 2(a). Subplot (b): The parameter values 
are as in subplot (a) with = 0.2, and for M = 5.3 (solid line), 5.4 (dashed line) and 5.5 (dotted line). Subplot (c): The 
parameter values are as in Fig. 2(c) with fii = O.f, and for M = 5.6 (solid line), 5.7 (dashed line) and 5.8 (dotted line). Subplot 
(d): The solid, dashed and dotted lines represent respectively for fid = 0.2, 0.4 and 0.6. The parameter values are Hi = 0.2, 
A'l = 5.4 and others as in Fig. 2(a). The points where V crosses the c/i-axis represent the amplitude of the SWs and the widths 
can be obtained as Iciml/IVJnini. 
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FIG. 5. (Color online) Numerical soliton solution (compressive) of Eq. (|13p for the parameters as indicated in the figure. The 
other parameters are as in Fig. 4(b). 
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FIG. 6. (Color online) Numerical soliton solution (rarefactive) of Eq. (|13|l for M — 5.2 (dashed line) and AI = 5.3 (dotted 
line). The other parameters are as in Fig. 4(c). 
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